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We study decays of D°, D + , and D~f mesons into two pseudoscalar mesons by expressing the 
decay amplitudes in terms of topological amplitudes. Including consistently SU(3 )f breaking to 
linear order, we show how the topological-amplitude decomposition can be mapped onto the stan¬ 
dard expansion using reduced amplitudes characterized by SU(3) representations. The tree and 
annihilation amplitudes can be calculated in factorization up to corrections which are quadratic in 
the color-counting parameter 1/N C . We find new sum rules connecting D + —> KsK + , Df —> Ksie + 
and D + —» K + n°, which test the quality of the 1/N C expansion. Subsequently, we determine the 
topological amplitudes in a global fit to the data, taking the statistical correlations among the var¬ 
ious measurements into account. We carry out likelihood ratio tests in order to quantify the role of 
specific topological contributions. While the SU(3 )f limit is excluded with a significance of more 
than five standard deviations, a good fit (with Ay 2 < 1) can be obtained with less than 28% of 
SU(3 )f breaking in the decay amplitudes. The magnitude of the penguin amplitude Pbreak, which 
probes the Glashow Iliopoulos Maiani (GIM) mechanism, is consistent with zero; the hypothesis 
Pbreak = 0 is rejected with a significance of just 0.7a. We obtain the Standard-Model correlation 
between B(D° —> Kltt 0 ) and B(D° —> A's7r°), which probes doubly Cabibbo-suppressed amplitudes, 
and find that B(D U —> Kltv 0 ) < B(D° — > Ksn°) holds with a significance of more than 4<r. We 
finally predict B{Df K l K+) = 0.012±°™® at 3cr CL. 


I. INTRODUCTION 


While there is a plethora of experimental information 
on hadronic charm decays, no theoretical method for dy¬ 
namical, QCD-based predictions for the corresponding 
decay amplitudes is known. The best theoretical ap¬ 
proach uses the approximate SU(3 )f symmetry of the 
QCD Lagrangian to relate the amplitudes of different 
decay modes to each other. If one assumes this sym¬ 
metry to be exact, one can express the amplitudes of 
all measured decay modes in terms of a smaller number 
of parameters, which are the reduced amplitudes char¬ 
acterized by SU(3 )f quantum numbers. Then one can 
predict the less precisely measured branching fraction 
on the basis of exact SU(3 )f or assess the validity of 
this assumption from the overall quality of the fit [Ml- 
SU(3)f is broken, because the masses m u ^d, s of the three 
lightest quarks are not equal. Comparing the differences 
among these masses with a typical hadronic scale one 
estimates SU(3 )f breaking to be around 30%. In prac¬ 
tice the quality of SU(3 )f symmetry can be much better 
(e.g. in heavy-hadron spectroscopy) or much worse (e.g. 
in heavy-quark fragmentation) and should be critically 
assessed for each system to which it is applied. Linear 
(i.e. first-order) SU(3 )f breaking can be rigorously in¬ 
cluded into the parameterization of the amplitudes, at 
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the expense of a larger number of reduced amplitudes. 
In the case of D —> PP', where D = D °, D + or Df 
and P, P' represent pseudoscalar mesons, such studies 
have been performed in Refs. [T5H2H]. (Remarkably, one 
can find relations between amplitudes which even hold 
to first order in SU(3 )f breaking (25].) Since there are 
fewer D —>■ PP 1 branching fractions than real parame¬ 
ters, there is a multi-dimensional space of solutions (all 
giving a perfect y 2 ) for the latter. Many of these solu¬ 
tions involve reduced SU(3)F-breaking amplitudes whose 
sizes are indeed of order 30% or less than the SU(3 )f- 
leading ones, giving evidence (but no proof) that the 
SU(3 )f expansion works. The redundancy associated 
with the multi-dimensional space of solutions poses a 
challenge for the numerical method to find the best-fit 
solutions because of the many flat directions in the space 
of reduced amplitudes. 

An alternative way to parameterize decay amplitudes 
involves topological amplitudes which are characterized 
by the flavor flow in the decays [51 IB HTTlUMTTUroi 12511501 - 
1521 . The building blocks of this approach are shown in 
Tab. [T] and Fig. [lj The topological amplitudes permit an 
easy and intuitive implementation of SU(3 )f relations. 
They further have the merit that they categorize the de¬ 
cays by dynamical criteria (i.e. whether the valence quark 
takes part in the weak interaction and which meson picks 
it up) and permit the combination of SU(3)i? methods 
with other calculational methods. In this paper we take 
a first step in this direction and apply the 1/N C expan¬ 
sion (first applied to D decays in Ref. [33]) to the tree 
(T) and annihilation (A) amplitudes of Tab. [lj (N c = 3 
is the number of colors.) T and A each factorize into 
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the product of a form factor and a decay constant up 
to corrections of order 1 /N 2 . We further include linear 
SU(3)j? breaking in the topological-amplitude decompo¬ 
sition, similar to the study of B decays in Ref. [T5]. For 
fixed values of T and A (obtained by adding a chosen 
1 /N 2 deviation to the factorized expressions) the num¬ 
ber of fitted complex topological amplitudes is reduced 
from 17 to nine, so that the problem of flat directions is 
substantially alleviated. 

The purpose of this paper is a systematic determi¬ 
nation of the topological amplitudes including linear 
SU(3) f breaking from a global fit to 16 D —> PP' branch¬ 
ing fractions and the measured strong-phase difference 
5k+tt~- For each topological amplitude we quantify the 
amount of SU(3)^ breaking with statistical likelihood- 
ratio tests using the statistical package rayFitter [53] . 
The latter is especially convenient in order to include 
nonlinear constraints in a frequentist analysis using the 
SLSQP algorithm implemented in SciPy [35] [331 • As 
a novel feature our statistical analysis fully includes the 
statistical correlations between the different experimental 
inputs. The ranges of the topological amplitudes found 
by us are an important input for the prediction of CP 
asymmetries. However, the latter also involve quanti¬ 
ties which cannot be extracted from branching fractions 
(SU(3)ir triplet amplitudes), so that additional input is 
needed for this purpose. This is one reason why we do not 
include measurements of CP asymmetries in our fit in¬ 
put. The other reason is their sensitivity to new physics, 
whose quantification should be separated from the deter¬ 
mination of hadronic parameters as much as possible. In 
this paper we also do not consider decays into final states 
with r) or rf, which involve additional parameters. 

The paper is organized as follows: In Sec. [H] we present 
the parameterization of D decay amplitudes using topo¬ 
logical amplitudes. We discuss the inclusion of linear 
SU(3 )f breaking and the appearing parametric redun¬ 
dancies in the diagrammatic language. In Sec. m we 
combine the method with l/N c counting and define our 
measures of SU(3 )f breaking. In Sec. IV we present the 
result of our fit. Finally, we conclude. 


II. DIAGRAMMATIC PARAMETERIZATION 
OF CHARM DECAYS 

A. Notation 


We choose the following conventions for the 

meson 

states: 



\K + ) = |us) , 

II 

(i) 

\K~) = - \au) , 

CO 

II 

(2) 

^ + 

II 

VT 

7T°) = -J= ( \dd) - \uu)) 

, (3) 

1 

II 

1 

^gj. 

D°) = - | cu) , 

(4) 

|P+> = | cd) , 

£>+) = |cs) . 

(5) 


Here the “=” sign means that the flavor quantum num¬ 
bers of the meson state on the left-hand side equal 
those of the quark-antiquark state on the right-hand 
side. Tab. [T] shows the topological (flavor-flow) ampli¬ 
tudes. The cross denotes the W-boson exchange encoded 
in the AC = 1 Hamiltonian. We write the Cabibbo- 
favored (CF), singly Cabibbo-suppressed (SCS), and 
doubly Cabibbo suppressed (DCS) decays as 


A GF (d) = v; s v u d A(d) = v c * s v ud <Ti , 

(6) 


A scs (d) = X ad A(d) = X sd Y, , 

(7) 


A DCS (d) = V* d V us A(d) = V c * d V us Y, . 

i 

(8) 

Here, 

we defined 



A s d := (A s - A d )/2 := (V c * V ua - V c * d V ud )/2 , 

(9) 


where A»d ~ A s — — A d . 7i is a topological amplitude (see 
Tabs Jll HI) and cf is the corresponding coefficient from 
Tab. [HI and d = D -A PP' labels the decay mode. There 
is a CKM-suppressed part oc V* b V u b in SCS amplitudes 
which can be safely neglected in all branching ratios. 

In the limit of unbroken SU(3 )f symmetry only the 
tree (T), annihilation (A), color-suppressed (C), and ex¬ 
change (E) amplitudes are needed to parameterize all 
D -A PP' decays. While the penguin amplitude P s ,d,b 
(labeled with the quark flavor running in the loop) is 
also non-vanishing in unbroken SU(3 )f, it only appears 
in the combination 


Pbreak = Ps ~ Pd, (10) 

where we have adopted the notation of Ref. [271. T,A,C, 
and E are commonly fitted together with the penguin 
amplitude Pbreak, which vanishes in the SU(3)f limit [9j- 
HU H31 HU 123 [37]. The normalization of the amplitudes 
is such that 


B(D -a PiP 2 ) = \A X {D -a PiP 2 )| 2 x V(D,P 1 ,P 2 ) , 


P(P,Pi,P 2 ) = t d x 


1 

16irm 3 D 



mp 2 ) 2 )(m 2 D 


( 11 ) 


(mpi +mp 2 ) 2 ). 

( 12 ) 


with X = CF, SCS, DCS. In the following, we will only 
make use of the notation A(d) without superscript, see 
Eqs. (§]§. 


B. SU(3)j?-breaking 

Any perturbative treatment starts with a subdivision 
of the Hamiltonian H = Hq + Hi into a piece Hq treated 
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TABLE I: SU(3)F-limit topological amplitudes. 


without approximation and the perturbation H±. The 
S-matrix element of the transition i f triggered by H\ 
is 

{f\Te~^ dixH ^ x) \i) . (13) 


In our case Hq is the QCD Hamiltonian with m u and 
m s set equal to m d - Hi consists of the weak |AC| = 1 
Hamiltonian Hw and the SU(3)F-breaking Hamiltonian 


Hvjj^f = ( m s - m d )ss, 


(14) 


where isospin breaking is neglected. With our choice of 
H 0 the asymptotic states i, f are eigenstates of H 0 which 
are D + or D° mesons or two-pion states. To first order in 
Hw and zeroth and first order in the transition 

amplitude in Eq. (131 becomes 


(/I - * J d 4 xH w (x) \i) + 

(f\ - 7 , JJ d^xd^THw^Hc^sjriy) | i) . (15) 

The second piece accounts for the differences of ampli¬ 
tudes involving a Df in the initial state or one or two 
kaons in the final state from their unflavored counter¬ 
parts. The Feynman rule of H^^j f is an ss vertex which 
we denote by a cross on the s-quark line. This approach 
is essentially identical to the one of Ref. |T9], where B 
decays have been considered. also leads to rj- 

rj' mixing. Using an 77 — 77 ' mixing angle in our diagram¬ 
matic method may lead to a double-counting of SU(3 )f- 
breaking effects and we do not consider final states with 


? 7 ^’s in the final state in this paper. The correspond¬ 
ing topological amplitudes are collected in Tab. [TIJ We 
combine our topological amplitudes into a vector 


P 


f f rri(l) rri(l) A a(1) a(1) <( 1 ) 

1 1 ) I | } JL 2 J J. 3 . /lj /IJ . X~L 2 5 713 ^ 

C,C^\C^\C^\E,E^\E^\E^\p hTe ^) T . (16) 


Then we can write 


M p = A 


(17) 


with a 17 x 17 coefficient 
(A(D°^K+K-),...,A(D+ 
the decay amplitudes. 


matrix M and A = 
-* K°K + )) subsuming 
The i-th column of M contains 
shows A in the 


III 


the coefficients cf of Eqs. (| 6 || 8 |). Tab. 
first column and lists the elements of M as table entries. 
We remark that the only final state with two identical 
mesons is | 7 r° 7 r°). In D°(po) n°( pi)tt°( p 2 ) two ef¬ 

fects must be taken into account: first, each topological 
amplitude appears twice (with pi and P 2 interchanged, 
leading to a proper Bose-symmetrized state). Second, in 
the subsequent phase space integration one integrates the 
azimuthal angle over the interval [ 0 , 7 r] rather the usual 
[0,27t], because the two pions are indistinguishable. The 
resulting factor of 1/2 in the decay rate (compared to 
the other listed decay rates) is accommodated thro ugh 
a factor of l/y/2 on the amplitude level in Tab. Ill For 
example, the factor of l/y/2 multiplying E is the result of 
the mentioned factors of 2 and l/y/2 and two factors of 
1 / 1/2 stemming from the | 7 T°) state in Eq. |3|. Note that 
it would be unwise to define the SU(3 )f limit from some 
average of s and d diagrams, since with this choice the 
asymptotic states constructed from H 0 would not corre¬ 
spond to physical mesons. Furthermore, there would be 


far fewer zeros among the coefficients in Tab. Ill which 
would further complicate the analysis. 

There is one more SU(3)F-breaking topological am¬ 
plitude, the penguin annihilation amplitude P^break = 
PA S — PA d depicted in Fig. [T] While the dynamics de¬ 
scribed by this amplitude is different from the ones dis¬ 
cussed so far, PAbreak enters the decay amplitudes in 
such a way that it can be absorbed into other ampli¬ 
tudes. Thus it is a redundant fit parameter, as explained 
in the following section. 


C. Redundancies 

The relationship between physical and topological am¬ 
plitudes is not one-to-one. If no other dynamical informa¬ 
tion on the latter is used, the determination of p from A 
in Eq. ( fl7| yields an infinite set of solutions describing the 
data equally well. A priori this feature renders fitted nu¬ 
merical values of T,..., Pbreak meaningless and obscures 
the comparison of different analyses in the literature. 
There are two ways to address this problem: one can 
simply remove redundant parameters and quote numbers 





























TABLE II: SU(3)F-breaking topological amplitudes corresponding to the amplitudes in Tab. Ill The Feynman rule for H^sj^gy F 
is the cross placed on an s line. 
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Decay ampl. A(d) 

T 

pW 

'T’(l) 
1 2 

rn(X) 

1 3 

A 


4 1} 

4 1J 

c 

a p) 


dP 

E 

4 L) 

e 

E ( 3 1] 

-Pbreak 

SCS 

A(D° -»• K+K~) 

1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

0 

1 

A(D° — > Tr + n~) 

-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-1 

0 

0 

0 

1 

A{D° ->• K°K°) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-1 

-1 

1 

0 

A(D° ->• 7r°7r°) 

0 

0 

0 

0 

0 

0 

0 

0 

1 

V2 

0 

0 

0 

1 

V2 

0 

0 

0 

1 

n /2 

A(D + ->• 7r°7T + ) 

1 

V2 

0 

0 

0 

0 

0 

0 

0 

1 

V2 

0 

0 

0 

0 

0 

0 

0 

0 

A(D+ -> K°K+) 

1 

1 

1 

0 

-1 

0 

0 

-1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

A{D+ -»• K°tv+) 

-1 

0 

0 

-1 

1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

A{Df -> K+ 7T°) 

0 

0 

0 

0 

1 

V2 

1 

\/2 

1 

\/2 

0 

1 

V2 

0 

0 

1 

\/2 

0 

0 

0 

0 

1 

v/2 


CF 


A{D° -)■ K~ 7T+) 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

0 

0 

0 

A{D° -»• K°n°) 

0 

0 

0 

0 

0 

0 

0 

0 

1 

\/2 

1 

V2 

0 

0 

1 

a/2 

1 

a /2 

0 

0 

0 

A(D+ ->• A'°7r+) 

1 

1 

0 

0 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

0 

A(D+ -> K°K+) 

0 

0 

0 

0 

1 

1 

0 

1 

1 

1 

0 

1 

0 

0 

0 

0 

0 


DCS 


A{D° -> K+ 7T-) 

1 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

1 

0 

0 

A(D° -> K° 7T°) 

0 

0 

0 

0 

0 

0 

0 

0 

1 

v /2 

0 

1 

a /2 

0 

1 

a /2 

0 

1 

\/2 

0 

0 

A(D+ K°n+) 

0 

0 

0 

0 

1 

0 

1 

0 

1 

0 

1 

0 

0 

0 

0 

0 

0 

A(D + K + n°) 

1 

a /2 

0 

1 

a /2 

0 

1 

a /2 

0 

1 

a /2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

A(D+ -» K°K+) 

1 

0 

1 

1 

0 

0 

0 

0 

1 

0 

1 

1 

0 

0 

0 

0 

0 


TABLE III: The coefficients of the decomposition of the physical amplitudes (including SU(3 )_f breaking) in terms of the 
topological amplitudes as in Eqs. (16] |8[) . The table entries are the elements of the coefficient matrix M in Eq. (171. 



FIG. 1: Penguin annihilation diagram. 


for the linear combinations of the topological amplitudes 
which are in one-to-one correspondence with the physical 
ones. Or one can use further theoretical (and experimen¬ 
tal) input to constrain the topological amplitudes. We 
determine redundancies among T ,..., Pbreak in this sec¬ 
tion and relegate the second approach to Sec. |III| 

The first redundancy is related to PAbreak of Fig. |T] 
which appears in SCS decays with the coefficients in 
Tab. |IV| The listed column of coefficients is linearly de¬ 
pendent on the four columns of coefficients of E, 3 
in Tab. |III| I.e. we can absorb P^break into the exchange 
amplitudes by redefining E = E — PAbreak, P-!/ 3 = 

p /2 3 + P^4break- In Tab. Ill this redefinition is implicitly 
already performed, so that PAbreak is not shown there 


Decay d 

P -/^break 

SCS 

D° -> K+K~ 

1 

D°—>n + n~ 

1 

D° -»• K°K° 

-1 

D° 7r°7r° 

1 

a /2 

D + — > 7T°7r + 

0 

D+ -> K°K+ 

0 

D+ K° 7T+ 

0 

D+ -> K+-k° 

0 


TABLE IV: The coefficients of the topological amp litude 
PA break which is absorbed into E, 3 in Tab. Ill as ex¬ 
plained in Sec. |II C| 


anymore. The physical meaning of E, E changes ac¬ 
cordingly, to be read as E, E[ 11 with the penguin anni¬ 
hilation mechanism included. However, |PAb re ak| is ex¬ 
pected to be negligibly small: The corresponding Wilson 
coefficient in H\\r is small and the momentum flowing 
through the penguin loop is large (of order of the D° 
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p 

ni 

n 2 

n 3 

n 4 

n 5 

ne 

p - Pbrkni 

p — An 4 

P ~ 4 U n 5 

T 

1 

1 

0 

1 

0 

0 

T — Pbi-k 

T - A 

T 

T w 

-1 

-1 

1 

0 

0 

0 

T[ 1] + P brk 

T w 

T w 

'T'(l) 

1 2 

-1 

0 

-1 

0 

0 

0 

P 2 (1) + Pbrk 

7^(1) 

1 2 

7■>(!) 

1 2 

7^(1) 

1 3 

0 

0 

0 

0 

0 

-1 

'T’(l) 

1 3 

n -t(l) 

1 3 

rpiX) 

-*3 

A 

0 

0 

0 

1 

-1 

0 

A 

0 

A + A« 

4 1} 

0 

0 

1 

0 

0 

-1 

A« 

A« 

A« 

4° 

0 

1 

-1 

0 

1 

0 

4 (1) 

A« 

4 (1) - 4 (1) 

^2 ^3 

A (1) 

0 

0 

0 

0 

1 

0 

A (1) 

^3 

4 1} 

0 

c 

-1 

-1 

0 

-1 

0 

0 

C + -Pbrk 

C + A 

c 

cP 

1 

1 

-1 

0 

0 

0 

cf ’ - Pbrk 

cp 

cf> 

cP 

1 

0 

1 

0 

0 

0 

- Pbrk 

cP 

dP 

Pi 15 

0 

0 

0 

0 

0 

1 


dP 

cP 

E 

0 

-1 

0 

-1 

0 

0 

E 

P +A 

E 

e\ 1] 

0 

1 

-1 

0 

0 

0 

eP 

P) 1} 

P( 1} 

eP 

0 

0 

1 

0 

0 

0 

pP 

P« 

eP 

EP 

0 

1 

0 

0 

0 

0 


E ^ 

^3 

eP 

Pbrk 

1 

0 

0 

0 

0 

0 

0 

Pbrk 

Pbrk 


TABLE V: The parameter vector p as defined in Eq. (161, 
vectors rq spanning the kernel of the coefficient matrix M 
and several redefined parameter vectors, see 


in Tab. 

III 

Eqs. (18 

20 


mass) so that the Glashow-Iliopoulos-Maiani (GIM) )38 
suppression will be effective. 

Further redundancies are related to the fact that our 
coefficient matrix M in Eq. (JTt| does not have maxi¬ 
mal rank. Considering first the SU(3)jr limit ignoring 
Tp, A- 1 ), C„-and eP one observes that the remain¬ 


ing matrix in Tab. |III| linking T, C, A and E to the 
physical amplitudes has only rank three. I.e. one of T, 
C, A and E is redundant. 

Redundancies of the diagrammatic approach in the 
SU(3 )f limit are also discussed in Ref. 5], comparing to 
the SU(3 )f parametrization in Ref. [3- The correspond¬ 
ing matching for B decays is done in Ref. [TO]. Note that 
the redundancies change when taking final states into 
account mna, leading to more parameters but also ad¬ 
ditional sum rules [251 . 

In clud ing SU(3)i? breaking, the 17x17 matrix M in 
Tab. Ill has rank 11. Consequently p in Eq. (16) con¬ 
tains six redundant complex parameters. The remaining 
parametric redundancy contained in Tab. Ill can be sys¬ 
tematically found and removed as follows. It is encoded 
in the six-dimensional kernel of the coefficient matrix. 
The 17-dimensional basis vectors of the kernel are given 
in columns 2 to 7 of Tab. [V] If we redefine p in Eq. (16) 
as 


p new = p + ^2 c.i n, : , a e C , (18) 

i 

this will not change M p in Eq. © , i.e. the n, define 


the “flat directions” in parameter space which correspond 
to the same One can remove this redundancy by 
redefining the topological amplitudes and choosing 11 of 
them as new independent parameters. For example we 
can set 


p nGW = P - ^breakni (19) 

which gives the result in the first column after the double 
line in Tab. [V] Subsequently, we can redefine the param¬ 
eters in order to eliminate Pbreak- In order to remove all 
redundancies in one step one can choose 

p new = p Pbreakni - EP U 2 ~ ' rl 3 “ 

An 4 — Ag 1 ^(n 4 + n 5 ) — C^ne (20) 

and then perform redefinitions of the other parameters 
in order to remove 


P hre ^,E^ ,E?\ A, A£> ,C, 


(i) 


( 21 ) 


from the parameterization. Note the special form of n 4 
which encodes the redundancy present in the SU(3)ir 
limit. n 4 forces us to eliminate one of T, C, E, or A, 
while the other five eliminations involve SU (3) ^-breaking 
amplitudes (e.g. those in Eq. ©)■ The elimination of A 
only is also shown in Tab.[V] Additionally, from n ,5 we see 
that the coefficient vector of Ag 1 is linearly dependent on 
the other annihilation coefficient vectors. Consequently, 
Ag 1 " 1 can be absorbed by redefining annihilation ampli¬ 
tudes only, as shown also explicitly in Tab. [V] 

Note further that the n; arc linearly independent also 
when removing all but the first six elements. This means 
it is not possible to perform redefinitions without touch¬ 
ing the tree or annihilation diagrams. Equivalently, the 
submatrix obtained by removing tree and annihilation 
diagrams from Tab. |IIl| has rank nine, which in this case 
equals the number of remaining parameters, i.e. the lower 
nine components of p. This observation guides us to the 
approach of Sec. |III| calculating tree and annihilation 
amplitudes will also remove the redundancies. 

After absorbing some topological amplitudes (e.g. 
those in Eq. (21)) into others the new amplitudes have 
lost their original meaning in terms of QCD dynamics. 
An important question in charm physics is the level of 
GIM cancellation between an s and d loop. In this pa¬ 
per we encounter Pbreak as a quantity probing the GIM 
mechanism. A naive quark-level calculation involves a 
suppression factor of m 2 s /m 2 and renders Pbreak negligi¬ 
bly small. Thus any information on the actual size of 
I Pbreak | may give insight into a possible non-perturbative 
enhancement of GIM-suppressed amplitudes. However, 
as shown above and exemplified in Tab. [V] the fit to 
topological amplitudes alone cannot give this informa¬ 
tion, because Pbreak cannot be separated from the other 
parameters fitted from the data. 

As we have seen above, the calculation of the kernel 
gives a method to remove redundant parameters. In the 
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same way, the cokernel of M gives us information on “re¬ 
dundant” amplitudes, i.e. six sum rules fulfilled by the 
latter, all of which were found in Ref. [25]. In other 
words: if one did a Gaussian elimination to determine p 
from Eq. one would end up with a 6x17 block of 
zeros in the transformed coefficient matrix M and linear 
combinations of physical amplitudes in the correspond¬ 
ing six entries of A. These linear combinations vanish 
by the SU(3) F sum rules of Ref. [55]. Thus the discussed 
redundancies are not the consequence of missing exper¬ 
imental information but of the symmetry relations un¬ 
derlying these sum rules. It is instructive to rederive 
these sum rules with our diagrammatic method, which is 
particularly straightforward and intuitive. We do this in 
Appendix [D] 

We checked that after the removal of all redundan¬ 
cies, the diagrammatic parameterization and the com¬ 
mon expansion in terms of SU(3) F representations can 
be mapped onto each other, i.e. one can calculate one set 
of parameters when given the other one. The mapping 
can be obtained explicitly by inverting either the reduced 
coefficient matrix M or its counterpart in the SU(3 )f 
method. Note that in the SU(3) F parameterization un¬ 
physical degrees of freedom are present in the very same 
way. Analogously, it is possible to redefine SU(3 )f ma¬ 
trix elements in order to obtain a physical basis 123] . In 
Appendix Blwe give the inverse of the SU(3 )f coefficient 
matrix of 23] and show the result of the extraction of the 
corresponding SU(3) F matrix elements for an example fit 
point of our diagrammatic analysis. So far our discussion 
of redundancies has assumed that the amplitudes in A 
are known. In practice, there is no information on most 
of their complex phases (and not all of them are physi¬ 
cal). This feature introduces additional flat directions in 
the space of our fit parameters and is equally present in 
the SU(3) F method. 

The discussion above has made clear that the 
topological-amplitude method is complete in the sense 
that it contains the full information contained in an 
SU(3 )f analysis including SU(3) F breaking to linear or¬ 
der. It is also worthwhile to study this question from 
the viewpoint of QCD dynamics: are there any dynam¬ 
ical mechanisms which cannot be mapped onto topolog¬ 
ical amplitudes? As a first topic we discuss final-state 
rescattering, i.e. decays D —»/'—>• / passing through an 
on-shell intermediate state /'. The flavor flow for such a 
rescattering process is always a deformation of a diagram 
in Tab. |T] or Fig. [l] and is therefore included in the cor¬ 
responding topological amplitude. Rescattering effects 
cannot be isolated from the “direct” D f decay, be¬ 
cause the dispersive part of A(D —► /' —> /) cannot 
be separated from that of A{D —> f) in a meaningful 
way. (Neglecting CP violation we can choose phase con¬ 
ventions such that the dispersive and absorptive parts of 
some amplitude equal its real and imaginary parts, re¬ 
spectively.) By the optical theorem the absorptive part 
of A(D —> f ) can be related to A(D —► /') and the 
f —> f scattering amplitude, with summation over all 



FIG. 2: Example for a SU(3) F -breaking diagram involving 
sea quarks which can be absorbed into Pbreak, see Tab. [IT] 


intermediate states f. This feature holds true for the 
topological amplitudes as well. The imaginary parts of 
the topological amplitudes found in our fit in Sec. [TV] 
are therefore a measure of the size of rescattering. The 
second topic of QCD dynamics addresses the proper de¬ 
scription of meson states. The state of e.g. an energetic 
kaon can be expanded as 

|/\°) = | ds) + \dsg) + | dsqq) + ... , (22) 


where the notation implicitly contains convolution inte¬ 
grals over the kaon momentum fraction carried by the 
indicated partons. Our graphical description of the topo¬ 
logical amplitudes only catches the first term in Eq. (221. 
The higher Fock states | dsg ), | dsuu ),... are suppressed 
with powers of the kaon energy, but in view of the small 
energy release in D decays this suppression is unlikely 
to be realized numerically. We may wonder whether the 
contributions with additional qq pairs in Eq. (22) will 
require the introduction of further amplitude topologies, 
with extra quark lines connected with “sea” quarks in 
the mesons. An example is shown in Fig. [2] However, it 
is easy to see that such diagrams are always obtained by 
forking a quark line of one of the topological amplitudes 
considered so far. For instance, the diagram in Fig. [2] is 
contained in Pbreak- 


III. THEORETICAL INPUT ON 
DIAGRAMMATIC SU(3) F BREAKING 


The great advantage of the flavor-flow parameteriza¬ 
tion over the plain SU (3) f approach is the opportunity to 
add dynamical input to constrain individual topologies. 
We use two different such inputs which are presented be¬ 
low. 


A. 1 /A c counting 

The 1/N C expansion [39: has first been applied to 
charm physics in Ref. [33] ■ We will apply 1/N C counting 
to the tree and annihilation topologies, which are leading 












Decay ampl. A(d) 

St — Sa 

C = C + S A 




E = E + Sa 

e[ 1] 

eP 


-Pbreak 

SCS 

A(D° -a K+K~) 

1 

0 

0 

0 

0 

1 

1 

1 

0 

1 

A(D° —> TT+TT”) 

-1 

0 

0 

0 

0 

-1 

0 

0 

0 

1 

A(D° -> K°K°) 

0 

0 

0 

0 

0 

0 

-1 

-1 

1 

0 

A(D° -A 7r°7r°) 

0 

1 

y/2 

0 

0 

0 

1 

x/2 

0 

0 

0 

1 

v/2 

A(D + -A 7T°7r + ) 

1 

x/2 

1 

0 

0 

0 

0 

0 

0 

0 

0 

A(D+ -A K°K+) 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

A(D+ -a K 0 7t+) 

-1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

A(D+ -A K+n°) 

0 

1 

x/2 

0 

0 

1 

x/2 

0 

0 

0 

0 

1 

x/2 


CF 


A(D° -A I<~ 7T+) 

1 

0 

0 

0 

0 

1 

1 

0 

0 

0 

A(D° -A K°n°) 

0 

1 

v /2 

1 

\/2 

0 

0 

1 

v /2 

1 

v /2 

0 

0 

0 

A(D+ -A K°n+) 

1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

A(D+ -a K°K+) 

0 

1 

1 

0 

1 

0 

0 

0 

0 

0 


DCS 


A(D° -a K+ 7 T-) 
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0 

0 

0 

0 

1 

0 

1 

0 

0 

A(D° -a K°tv°) 

0 

1 

v /2 

0 

1 

v /2 

0 

1 

v /2 

0 

1 

v /2 

0 

0 

A(D+ -A K°n+) 

0 

1 

0 

1 

0 

0 

0 

0 

0 

0 

A(D+ -A K+n°) 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

A(D+ -A K°K+) 

1 

1 

0 

1 

1 

0 

0 

0 

0 

0 


TABLE VI: Coefficients of the paramet ers ( St — Sa , . . • , TWeak) for the amplitudes .4(d), which are obtained from .4(d) by 
subtracting the factorized part, see Eq. (27l. The table entries are the elements of the coefficient matrix M in Eq. (34). 


in l/iV c . Here we exemplify the method for T: 


tfac 


T =T 
Gf 

V2 


nfac 


' + St, 

(23) 

aifn(m 2 D - ml)F^(ml ). 

(24) 


Here ai = C 2 +C\/N c = 1.06 in terms of the usual Wilson 
coefficients C \ t 2 of Hw and the quoted value corresponds 
to next-to-leading order in the NDR scheme at a scale of 
1.5 GeV. It is important to note that the color exchange 
between the two quark lines in the T diagram in Tab. [I] is 
penalized by two powers of 1/IV C . We parameterize thi s 
1 /N% correction by the complex parameter St in Eq. ( p3| ). 
Also the renormalization scale and scheme dep end ences 
of ai are suppressed by 1 /N 2 . By using Eq. (23) and 
the equivalent formulae for the other tree amplitudes 
T + ... we trade four parameters for a single pa¬ 

rameter St with |<5T/T fac | < 0.15. SU(3 )f breaking in 
this small parameter is neglected, because it is smaller 
than the neglected second-order SU(3)^-breaking effects. 
F^ >n (m'^) entering T fac is measured in semileptonic D de¬ 
cays, therefore the 1 /N c method uses additional exper¬ 
imental input, too. Also the A amplitudes factorize up 
to corrections of order 1 / N 2 . The factorization formulae 
for all tree and annihilation amplitudes can be found in 


Appendices |C 1| and |C 2[ respectively. In analogy to St 
we define the complex parameter 


5 A = A - A 


fac 


(25) 


for the 0(l/N 2 ) corrections. A fa c dep ends on the form 


factor Ej f7r (m| J(s) ), see Appendix 


C 2 


for details. 


E, C , and Tbreak are formally suppressed by one power 
of 1 /N c with respect to T. However, E and C are en¬ 
hanced by short-distance QCD effects residing in the 
Wilson coefficients: we write Hw oc C\Q\ + C 2 Q 2 = 
(Ci + C 2 /N c )Q 1 + 2 C 2 Q& with the octetxoctet operator 
Q s = u'y ll T a cq l j ,1 T a q and note that (PP'\ Q 8 \D) is 1 /N c 
suppressed. However, the Wilson coefficient 2 C 2 = 2.4 
almost exactly offsets the l/iV c suppression, so that E 
and C can be almost as large as T. We therefore do not 
place a numerical constraint on \E\, \C\, or |Pbreak| in 
our fit but rather keep them general. 

With the added 1/N C input the diagrammatic anal¬ 
ysis becomes more constrained compared to the plain 
SU(3)j? approach. Factorization fixes the sizes of the 
tree and annihilation amplitudes within roughly ~ 15% 
of T fac , i.e. the size of the 1 /N 2 corrections. In the case 
of A the 1/N% corrections quantified by Sa include final- 
state rescattering effects HUSK which are not pro¬ 
portional to the decay constant fo which enters A fac . 
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We therefore do not normalize 5 a to A fac , but instead 
allow \5a to be as large as |5 t|- Factorization has also 
been used in Refs. [2B1 02], but only to estimate SU(3 )f 
breaking. We instead use it to constrain the overall sizes 
of T and A. Note that we treat T lac and A fac beyond lin¬ 
ear SU(3 )_f breaking, so that these factorized amplitudes 
violate the Grossman-Robinson SU(3)^ sum rules [25| . 

The parameters St and 5a replace the first eight entries 
of p in Eq. (16) as fit parameters. We use 


p ' = (s T ,5 A ,C,C^,C^\c, 


( 1 ) 

3 J 


E, P^P^P^Pbreak) 


(26) 


comprising 11 parameters in total. We next derive the 
equivalent of Eq. <0 for this new set of parameters. To 
this end we define 


A(d) = A(d)~A iac (d), 
A iac (d ) = T {ac (d) + A tac (d). 


(27) 

(28) 
to 


The 17 x 11 coefficient matrix linking p' 

A = (A{D 0 -»■ K+K~),...A(D+ -t K 0 K+)^j has 
only rank 10. This has two implications: First, there is 
still a redundant parameter. Second, there is a new sum 
rule among the physical amplitudes. Addressing the first 
point, the kernel has the 11-dimensional basis vector 


n= (-1,-1,1,0,0,0,1,0,0,0,0)' 


(29) 


where the order of the entries is the same as in Eq. (261. 
The redefinition 


P" = P' 


5 ,in 


(30) 


with n as in Eq. (29) absorbs 5a into C, E and 5t' 
setting 


C = C 
E = E 


5 a , 
5 a , 


(31) 

(32) 


one observes that the physical amplitudes only depend 
on C, P, and 5t — 5a- Writing 


p = (St-6 a ,C,c£\c£\ 

r>(l) 77>(1) 7T’(l) 


Cl 


( 1 ) 


TA 7 A\ ± J JAX^J T7'\ ± J P 

5-^2 5-^3 5 -‘break 


(33) 


the desired equivalent of Eq. (17) reads 


Mp = A , (34) 

with the amplitudes of Eq. ( [27]) o n the RHS. The result¬ 
ing 17 x 10 coefficient matrix M with rank 10 is shown 
in Tab. |VlJ Addressing the second point, we find the new 
sum rule 

A(D+ K°K+) - A(D+ -» K°n + )- 

2V2 A(D+ ->• K + tt°) = 0, (35) 


from the cokernel of M. It tests the 1/N C counting and 
is violated by terms which are linear in SU(3 )f breaking, 
but suppressed by two powers of 1/N C . 

The 1 /N 2 corrections parametrised by 5t,a are varied 
in smaller ranges than the other fit parameters. If we 
consider them fixed, there remain nine unknown complex 


parameters in Eq. (261 and a corresponding coefficient 
matrix with rank nine, implying a new sum rule. We 
combine the new rule with the one in Eq. (|35|) as: 


A(D+ -► K°K+) - A(D+ -»■ K°ir + ) 

= 2 (5t-S a ), (36) 

A{D+ -»■ K+tt°) = 1=(6t- 5a) ■ (37) 


The amplitudes in Eq. (35) are related to those with 
Ks,l in the final state as 

A(D+ K s , l K+) = A(D + K°K +), (38) 

v 2 


A(Df K s ,l7t + ) = -j=A (D+ -»> A°7t + ) . 

The corresponding branching ratios read: 

B(D + -> K s , l K+) = \\ sd \ 2 V(D+,K°,K+)x 
| A faC (D+ -»■ K°K+) + ( S T - S A ) + Pbreakl 2 
B(D+ K s ,l7t + ) = \X sd \ 2 r(Dt,K°,7r + )x 
| A iac (Dt -> K° 7T+) - (S T - <5 a ) + P break | 2 , 
B(D+ K+ 7T°) = \V^V us \ 2 V(D+, K+,n°)x 
\A lac (D+ ->• K+ 7T°) + (S T - <5a)| 2 , 


(39) 


(40) 

(41) 

(42) 


with P(P>,Pi,P 2 ) as defined in Eq. (12). Eqs. (40-42) 
permit one to probe our combined SU(3 ).f and 1 /N c ex¬ 
pansion quantitatively, since a too large value of |<5-r — 5 a \ 
extracted from Eqs. (40-42| would falsify the method. 
Furthermore, the size of |Pbreakl gives insight into an 
important issue of QCD dynamics, the size of the GIM 
suppression in the difference between strange and down 
loops. 

In Tab. [YTI 


we show example fits to the branching 


ratios Eqs. (40) (42) only, testing the dependence of the 
fit result on the 1 /N£ corrections and the broken penguin. 

In the first place, we illustrate that the data can easily 
be accommodated for realistic values of St, 5a and P br eak 
(point 0). Taking these parameters out of the fit (point 1) 
results in a bad description of the data which is rejected 
at 7a. It is possible to describe the data with an enhanced 
broken penguin only (point 2) but also with P bre ak = 0 
and adjusting St and S A (point 3). A better kno wledge 

), see Appendix 


of the form factor 


’(r 


D. 




C 2 


is cru¬ 


cial in order to disentangle an enhanced penguin from 
1 /N 2 corrections. This could be provided by future high 
statistics measurements of r decays (HO [50]. 
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Example fit point with minimal y 2 : 

Point 0 

Point 1 

Point 2 

Point 3 

Exp. data 

Applied conditions: 

None 

-^break = $T = $ A — 0 

o 

II 

II 

h 

-^break ~ 0 

— 

|Pbreak/T faC | 

0.25 

0 

0.54 

0 

— 

arg(Pbreak) 

-1.95 

0 

2.21 

0 

— 

\5 T \/T lac 

0.11 

0 

0 

0.15 

— 

arg (<5 t ) 

3.07 

0 

0 

-2.74 

— 

\ 5 a \/ A { ^ 

0.09 

0 

0 

0.11 

— 

arg(<bi) 

0.67 

0 

0 

0.00 

— 

F q d ’ k (0)/F o d A0) 

0.96 

1.01 

0.95 

0.97 

— 

F 0 DK (0) 

0.74 

0.72 

0.74 

0.74 

— 

Fo >7r (0) 

0.64 

0.64 

0.64 

0.64 

— 

IPf-(m 2 D(s) )| 

2.39 

1.99 

4.50 

1.62 

— 

arg 

1.71 

-1.29 

-1.15 

-2.36 

— 

T fa C /i 0 -6 Ge y 

2.52 

2.52 

2.52 

2.52 


T fac (T>+ A'°A' + )/10“ 6 GeV 

3.40 

3.34 

3.40 

3.40 

— 

T fac (P+ A'°7r+)/10" 6 GeV 

-2.53 

-2.68 

-2.51 

-2.57 

— 

T fac (A>+ ->■ A' + 7r°)/10" 6 GeV 

2.22 

2.22 

2.22 

2.22 

— 

A Cac (Df -r A'°7r + )/10 -6 GeV 

-0.18 + i 1.22 

0.28 -i 0.99 

0.94-*2.12 

-0.59-*0.59 

— 

A f ac(L>+ A' + 7r°)/10 -6 GeV 

0.10 -i 0.68 

-0.16 + *0.55 

-0.53 + *1.19 

0.33 + *0.33 

— 

B{D+ A s A + )/10“ 3 

2.83 

4.04 

2.85 

2.83 

2.83 + 0.16 [46] 

B(D+ ->■ Astt + )/10~ 3 

1.22 

1.22 

1.23 

1.22 

1.22 + 0.06 tglSS] 

B(D+ -a A + 7r 0 )/10~ 4 

1.83 

1.83 

1.72 

1.83 

1.83 + 0.26 [46] 

,x 2 

0.00 

63.93 

0.20 

0.00 

— 

V 

— 

5 

3 

2 

— 

Significance of rejection 

— 

7.0a 

0.03cr 

0.0(7 

— 


TABLE VII: Fits to the branching ratios B(D + —> I\sK + ), B(D+ —> Kg 7r + ), B(D + —> K + n°) only, without taking correlations 
and additional constraints on SU(3 )f breaking, see Sec. |IIIB| into account. The form factors are varied as described in Sec. IV 
and Appendices C 1 and C 2 The y 2 is the one taking into account the three given branching ratios and form factors only, v 


are the number of degrees of freedom compared to the fit scenario of point 0. The significance of rejection takes point 0 as null 
hypothesis. ^Our average. 


B. Measuring diagrammatic SU(3 )f breaking 


In order to describe SU(3) f breaking in the framework 
of the diagrammatic approach, we introduce the following 
measures in analogy to Ref. [27. We define 


5 


f,T 

X 


= ma Xd 


•Axjd) 

A(d) 


(43) 


where T = C, -E^Pbreak and Ax(d) is the part of the 
amplitude of decay d stemming from the corresponding 
SU(3)F-breaking parameter(s) only. A(d) denotes the 
full amplitude of decay d. The parameters defined in 
Eq. ( Jl3j ) give a measure for the maximal SU(3)F-breaking 
contribution to the full amplitude from each topology. 

A measure of the maximal SU (3) f breaking residing 
in any of the topologies C, E and Pbreak is therefore 


cV,topO 

Ox F = max d 


Sr -^x(d) 

A{d) 


(44) 


Note that the SU(3 )f breaking stemming from our 
calculation of the T and A topologies using factorization 
is not included in the definition Eq. (44). 


Furthermore, we quantify the relative SU(3 )f breaking 
of C and E topologies by the measures 


,cfVc 


a 


(i) 


c 


s: E i 1) / E 
°X 


E, 


(i) 


E 


(45) 


respectively. In the fit we always demand all the above 
measures to be < 50%. In S'p and <5^ topo we ignore 
B(D° —> K°K°) when taking the maximum, because this 
branching ratio vanishes in the SU(3 )f limit. Note that 
Eg 1 "* appears in the omitted channel D° —> K°K° only, 

therefore S'j^ = 0 is insensitive to the size of E^p ^ 0. 
Furthermore, in case an amplitude vanishes at some point 
in parameter space we also exclude it from the calculation 


of the maxima in Eqs. (43) and (44). 
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IV. FIT TO BRANCHING RATIO 
MEASUREMENTS 


In our global fit we use the available measured branch¬ 
ing fractions and the strong phase difference Sj^+n- and 
impose the theoretical constraints quoted in Sec. |III| 18 of 
the fit parameters are related to topological amplitudes: 


I. arg(C), 
I, arg(cf } ), 

|Pbreak/T faC |, 


|cyr fac 


\E/T { -% arg (E), 
\4\ arg (Pf } ), 

ai'g (/’break ), 


with i = 1,2,3 and T fac is calculated from Eq. (24). We 
normalize to T fac rather than T = T fac + St, because 
our fit is only sensitive to the combination St — 5a and 
therefore leaves St undetermined. These 18 quantities 
are supplied by four parameters measuring the 1/A/? cor¬ 
rections to the tree and annihilation diagrams: 


|<M/T fac , arg (S T ), \S A \/T iac , arg (S A ). 


In addition we need five parameters related to form fac¬ 
tors: 


F“- K {0)/F°*{ 0), 
\F 0 K ”(m 2 D )\, 


F 0 DK ( 0), ^(0), 

arg (F<F*(m 2 D )), 


and set F^' K (m 2 Ds ) = F^*(m 2 D ). Altogether these are 27 
real parameters, which are fitted to 16 measured branch¬ 
ing ratios and one strong phase. The experimental input 
values, including the respective correlations, are listed in 
Appendix [A] The number of parameters is larger than 
the number of observables. However, the 27 parame¬ 
ters are subject to 10 constraints on the maximal size of 
linear SU(3 )_f breaking, see Sec. IIIB and the bounds 
| Jr, a | < 0.15 T. At the global minimum we obtain 
X' 2 = 0.0, i.e., the parameterization and theoretical input 
is in perfect agreement with the data. Thus the data 
are both compatible with our chosen bound on SU(3 )f 
breaking (i.e. all measures defined in Sec. IIIB are smaller 


than 50%) and the six Grossman-Robinson SU(3 )f sum 
rules [291 . 

In order to study the relative importance of the topo¬ 
logical amplitudes for the description of the data, we per¬ 
form likelihood ratio tests. We look at several scenarios 
where some of the parameters of our fit are fixed. In order 
to keep the fit simple, we assume the validity of Wilks’ 
theorem ED, i-e-, we calculate the p-value according to 

[321161 


P= 1 - P v /2 (Ay 2 /2) 


(46) 


with the normalized lower incomplete Gamma function 
Pv /2 depending on the number v of relatively fixed pa¬ 
rameters compared to the full fit. For a general discussion 
of the assumptions underlying Eq. (461 see Ref. |34]. 


The results of our likelihood ratio tests are shown in 
Tab. |VIII| This table shows at which significance we can 
reject a certain hypothesis. For example, we can reject 
Pbreak = 0 at only ~ 0.7cr, implying that P br eak = 0 is 


well consistent with the data. However, the fit shows 
a clear need for SU(3 )_f breaking: the SU(3) p-limit fit 
with Pbreak = E = C^ = OVi is rejected at > 5a. 
Looking at the SU(3)j? breaking in specific topological 
amplitudes we find a slight tendency toward a stronger 
SU(3) f breaking in the color-suppressed tree than in the 
exchange diagrams. 

In Figs. !□ we show plots of the fit parameters, mea¬ 
sures of SU(3 )_f breaking and fit predictions for observ¬ 
ables. We see that in the multi-parameter space the best- 
fit solutions cover broad regions and typically several dis¬ 
connected best-fit regions exist. Considering that there 
are more parameters than fitted quantities the large de¬ 
generacy of the best-fit region is not surprising. It is 
moot to quote best-fit values for the parameters, because 
one can move in a wide valley with A\ 2 = 0. We sus¬ 
pect that the alternative approach of a Bayesian analysis 
would single out a small portion of this Ay 2 = 0 val¬ 
ley as a consequence of the Bayesian prior placed on the 
fit parameters and the central limit theorem of statis¬ 
tics. Therefore frequentist analyses like ours are more 
adequate to the problem. 


The phase of C significantly deviates from 0 and ±7r 
(see Fig. |3(b)| ), which points to large rescattering ef¬ 
fects. The fit results for |C'/T fac | and \E/T tac \, see Fig. [ 3 J 
show disconnected regions at 95% CL. C and E are sup¬ 
pressed by 1/N C but involve a large Wilson coefficient 
~ 2.4 as discussed in Sec. |III A| Thus only solutions with 
|C/T fac |, \E/T iac \ < 1 are consistent with 1/N C counting, 


which singles out one of the three regions in the |C/T fac |- 
\E/T tac \ plane. 

The needed maximum size of total SU(3 )f breaking 
on the amplitude is given by (5^ topo ~ 30% in agreement 
with Ref. [23], as can be read off Fig. |4(d)| Note that 
S'jtopo ag ag j//r n0 ^. measure the average but the 
maximal size of SU(3)^ breaking in one of the 17 decay 
channels except for D° —» K°K°, see Sec. IIIB Thus, 


these measures are very conservative and could in princi¬ 
ple be biased by a single channel. However, the need for 
SU(3) f breaking in individual parameters can be consid¬ 
erably smaller than 30%, e.g., |£’| 1 ^| ~ 0, IC^I ~ 0 is 
well allowed at la, see Figs. |5(a) and 5(c)| respectively. 
Also a | Pbrea k/T| bel ow 5 % already gives very good fits, 
see Fig. |6(a)| In Fig. |4(c)] we see that the same is consis¬ 
tently the case for S'^ k . 

From Figs. |4(a)| and |4(b)| we again see the slight ten¬ 
dency for larger SU(3 )f breaking in the color-suppressed 
tree topologies compared to the exchange diagrams. 

As illustrated in Tab. IVIIl discussed in Sec. lIII Al around 


Eqs. (|40j|42j) , the broken penguin Pbreak is correlated 
with the parameter St — S A quantifying 1 /N 2 corrections 
to factorizable amplitudes. This feature can be verified in 
Fig. |6(c) which shows this correlation. A vanishing pen¬ 
guin Pbreak ~ 0 is allowed at the price of 1/N C breaking 
corrections of order > 15%. Note again that this correla¬ 
tion heavily depends on the poorly measured form factor 
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Hypothesis 

Significance of rejection 

AX 2 

dof 

Rbreak 0 

0.7cr 

1.3 

2 

Pbreak = = Cf > = OV* 

> 5(7 

431.4 

14 

= 0 Vi 

3.0cr 

20.3 

6 

E = E < 1 1) * =0 Vi 

> 5(7 

156.4 

8 

cf 5 = OVi 

4.3(7 

31.6 

6 

C = cf) = ov* 

> 5(7 

267•10 3 

8 


TABLE VIII: Results of several likelihood ratio tests. Shown are the obtained y 2 , the relative degrees of freedom ( dof ) of the 
hypothesis compared to the null hypothesis, which is the full fit, and the significance at which the hypothesis can be rejected. 


F^™(m 2 D } ). Interestingly, the fit result for ( rn D (B) )> 
see Fig. [7] is not completely flat, showing its nontrivial 
influence on the branching ratios of charm decays. The 
branching ratio B(D + —> K + n°) depends on no topolog¬ 
ical parameters besides St and 6a ■ Its fit result, which 
is given in Fig. |8(a)| shows that our assumptions on the 
ranges for St and 5a are loose enough to accommodate 
the measured branching fraction. However, large fit re¬ 
sults for B(D + —> K + n°) are slightly disfavored. 


We may next ask whether we can use our fit out¬ 
put to predict individual branching fractions better than 
they ar e cur rently measured. Our general finding is as 
in Fig. 8(a) the fit output for the Ay 2 profiles essen¬ 
tially tracks the fit input. To find non-trivial predictions 
for future measurements we must study correlations be¬ 
tween at lea st two observables. A nice result is shown 
in Fig. 8(b) revealing the correlation of B(D° —> Kltt°) 
and b(fF Ksn°). In the SU(3 )f limit the branching 
ratios are strongly correlated through their parametric 
dependence 1 


B(D° ->■ K, s7T°) ~ \E - Cf + 2A 2 |E - Cf , (47) 
B{D° K l tt°) ~ | E- Cf - 2A 2 |E - C \ 2 , (48) 


findings further we define 


0 = B(D° -> K s n") - B(D° -f K L n°) 
1 ’ B(D° -»• K s n°) + £(L>° -»• K l tt°) 


(49) 


and quote the confidence intervals in the first row of 
Tab. |IX[ The ratio of the magnitudes of the DCS and 
CF amplitudes is listed in the third row of this ta¬ 
ble. Fig. [9] visualizes these confidence intervals and also 
shows the prediction of Refs. [521155] , which is the black 
dot corresponding to R(D°) = 2 tan 2 9c (where 6c is 
the Cabibbo angle). The result is quoted without un¬ 
certainty in these papers and Refs. [55H55] argue that 
corrections from SU(3) breaking to these relations are 
small. Refs. [54] [55] arrive at this conclusion by cal¬ 
culating the amplitudes in QCD factorization [sa na, 
which is a calculational method valid for values of m c 
much larger than the hadronic scale governing the in¬ 
frared structure of the decays. Our fit permits sizable 
corrections to R(D°) = 2 tan 2 6q from the SU(3) break¬ 
ing contributions, so that future measurements will give 
insight into the size of SU(3) breaking and the viability 
of QCD factorization in charm physics. 

Another test of doubly Cabibbo-suppressed contribu¬ 
tions involves the decays Df —>• K$,lK + . We study 


which implies B(D° —> K^n 0 ) < B(D° —► Kstt 0 ). 
This relation is a priori absent once SU(3)F-breaking ef¬ 
fects are included, because the latter can be larger than 
|E — C\. However, the global fit rejects this possibil¬ 
ity: in Fig. |8(b)| the region corresponding to 95% CL 
entirely satisfies B(D° —> Kl 7 t °) < B{D° -A- Ksn°). 
Performing a dedicated likelihood ratio test we find that 
B(D° -> K L n°) < B(D° -> K s n°) holds with a signifi¬ 
cance of more than 4er. Our fit excludes a large region of 
the B(D° —> Ksn°)-B(D 0 —> Kl 7t°) plane which is still 
allowed by the individual measurements. To quantify our 


R(Df) 


B(D+ ->■ K S K+) - B(£>+ -> I< l K+) 
B{D+ -»• K S K+) + B{D+ -> K l K+) 


and predict the not yet measured observables R(D+) and 
B(D+ —> KlK + ), see Tab. IX and Fig. [ 9 ] Agai: 


iin com¬ 


paring our result with the prediction in Ref. m we find 
much larger uncertainties. Thus also in Df —► Ks,lK + 
future data will test the accuracy of QCD factorization 
assumed in Ref. 1551. 


V. CONCLUSION 


1 In order to find the correct relative signs in Eqs. © and ( ] 4 8 [l one 

must define Ks,l correctly. Eqs. |Tj:i and (Dj) comply with | K ' ] ) = 

C\K°) = —CP \ K°) entailing \K S ) ~ (|A'°) - \K°))/y/2. We 

have checked our results by studying the full decay chain D° -4 

> 7r+7r — ]7r°, from which the K° sign conventions drop out. 


We have studied the decay amplitudes of D mesons 
into two pseudoscalar mesons with the topological- 
amplitude approach. To this end we have incorporated 
linear SU(3)jr breaking into the method and have shown 
that the topological amplitude method can be mapped 
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onto the standard decomposition of the decay ampli¬ 
tudes in terms of reduced amplitudes characterized by 
SU(3 )f representations. Unlike plain SU(3)^ analyses 
the topological-amplitude method permits the use of a 
1 /N c expansion to calculate the factorizable tree and an¬ 
nihilation amplitudes in terms of form factors and de¬ 
cay constants, up to corrections of order 1 /Nf.. This 
additional theoretical input has led us to a new sum 
rule between the branching fractions of D + — > K$K + , 
Df — > Ksir + and D + — > K + ir°. This sum rule corre¬ 
lates the non-factorizable 1 /V 2 terms with the penguin 
amplitude Pbreak • The latter quantity is of prime interest 
to understand the dynamics of flavour-changing neutral 
current transitions in the charm sector, because Pbreak is 
suppressed by the GIM mechanism and vanishes in the 
limit m s = rrid- 


We have then performed a global fit using all avail¬ 
able branching ratios and the experimental information 
on the strong phase difference Sj^+ n ~. In our analysis we 
have included the information on correlations between 
experimental errors. It is possible to find a perfect fit, 
with a large parameter region satisfying y 2 = 0. This 
means that current data comply with i) the Grossman- 
Robinson sum rules [29], ii) our chosen upper bound of 
50% on SU(3 )f breaking, and iii) our assumption that 
the 1/iV 2 corrections to the factorizable amplitudes are 
smaller than 15% of the factorized tree amplitude. The 
main phenomenological results of our paper are various 
likelihood ratio tests addressing the sizes of the topolog¬ 


ical amplitudes and their SU(3 )f breaking (Tab. VIII 
and Figs. 


_ j. Importantly, we find that there is no ev¬ 
idence for an enhanced broken penguin. The hypothesis 
Pbreak = 0 is rejected at below la only, i.e., insignifi¬ 
cantly. Improvements of B{D + —>■ K$K + ), B(Df — t 
A's , 7t + ), B(D+ —> K + tt°) and especially the form factor 
P ( f- 7 r (m| )( } ) could advance our knowledge of the GIM 
mechanism in charm by pinning down the proportions of 
broken penguin and 1/7 V 2 cor rections. The current sta¬ 
tus is summarized in Fig. 6 (c) While the SCS branching 
ratios B(D + —> K$K + ) and B(D+ —> A's 7 t + ) are known 
at a precision of < 6 % the relative uncertainty of the 
DCS branching ratio B(D + —> A' + 7t°) is about ~ 14% 
and leaves room for improvement. As the latter is the 


only charm decay into kaons and pions which depends 
on factorizable contributions only, it is very important 
to improve its measurement. With a simultaneously im¬ 
proved P ( f r 7 r (m| )(s) ) the branching ratio B(D + -± A + 7 t°) 
serves as a test of factorization in charm decays. 

We observe a slightly larger SU(3 )_f breaking in color- 
suppressed tree than in exchange diagrams. In no chan¬ 
nel more than ~ 30% SU(3 )_f breaking is needed to de¬ 
scribe the data (not considering P>° —> KgKs, which is 
forbidden in the SU(3 )f limit); this finding agrees with 
the plain SU(3 )_f analysis of Ref. [23,]. However, as a 
matter of principle one cannot decide whether or not the 
actual SU(3 )f breaking is larger than this. This can po¬ 
tentially only be achieved by future QCD calculations on 
the lattice [58] . In the data, there is no indication of this 
to be the case. 

With our topological-amplitude fit it is further pos¬ 
sible to make predictions for branching fractions which 
can be probed by future measurements. Despite our 
conservative ranges for the SU(3)^ breaking parame¬ 
ters, we find a correlation between B(D° —> K^n 0 ) and 
B(D° —)• As 7 t°) probing the doubly Cabibbo-suppressed 
contributions to these modes: Fig. [ 8 ] entails the predic¬ 
tion B(D° —> < B(D° —> A 57 T 1 ) at more than 4er. 
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Observable 

±1(7 

±2(7 

±3(7 

R(D°) 

R(D+) 

B(D+ -± K l K+) 

\A DCS {D°^K°n 0 )\ 


n OQ +0 04 

u.uy_ 00 2 

n 11 + 004 
0-012“ 
0.05tH? 

o.oslH 2 

n OQ + 0 - 07 
U.UJ-q 04 

n ii +0 06 

0-012“ 

0.05±°°3 

0.08«°6 

0 QQ+0.09 
U.UJ—005 

n 11 + 006 
u>11 -0.20 

0 012 +0 006 
W.U±Z,_Q 002 

0-05l°;°3 

0.08l°;° 4 7 

\ _A CF (D°-s.iC°7r°) 

A dcs (D+^K°K+) 

A cf (D+^K°K+) 


TABLE IX: Fit results for several observables probing doubly 
Cabibbo-suppressed amplitudes. The corresponding plots are 
shown in Fig. [9] 
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(c) (d) 




(e) (f) 


FIG. 3: SU(3)f limit topologies. In Figs, (e) and (f) the dashed (solid) line denotes the 68% (95%) C.L. contour. 
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FIG. 4: Ax' 2 profile of the parameters S'^' E ’ PhrBa - ] ‘ measuring SU(3)F-breaking in C, E, and Pbreak (a,b,c) and of S'x° po defined 
in Eq. (|44[) , which quantifies the overall SU(3)F-breaking (d). 
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icf’/ci 

(a) 


'X 

< 


i - 1 
i = 2 
* = 3 


Ul 


\l 




—7r —7r/2 0 7r/2 


arg(Cr 

(b) 




(c) 


(d) 


FIG. 5: SU(3)_F-breaking color-suppressed tree (a,b) and exchange (c,d) topologies. 





(a) (b) (c) 


FIG. 6: The broken penguin (a,b) and its correlation to parameters measuring the 1/IV C corrections (c). In Fig. (c) the dashed 
(solid) line denotes the 68% (95%) C.L. contour and the region to the right of the contours is allowed. 
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\Fo K% (m 2 D J 


FIG. 7: The form factor F,f*{m 2 D(g} ). 




(b) 


FIG. 8: Ax 2 profile of B{B+ ->• K+n°) (a) and correlation between B(D° -»■ K L n°) and B{D° -»■ TF S 7r°) (b). In Fig. (a) the 
red dashed (solid) line indicates the Ict (2ct) experimental error. In Fig. (b) the dashed (solid) lines are the 68% (95%) C.L. 
contours of our fit and the dark (light) gray shading denotes the 68% (95%) C.L. region of the measurements. Here, the solid 
red line corresponds to B(D° —¥ Kltt°) = B(D° —> Ksn 0 )- 
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-0.05 0.00 0.05 0.10 0.15 0.20 0.25 
R(D°) 



0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.10 

\A dcs (D° -> K°ir°)/A CF (D° ->• 


(a) Black: prediction of Refs. |52H55| . Red (below the black 
point): experimental error 


(b) Black: prediction of Refs. |52H55| 







1 -. .-1 

1 #1 



1 - 0 - 1 

1*1 


1*1 


R(Dt 


0.10 0.15 0.20 


0.00 0.02 0.04 0.06 0.08 0.10 0 . 12 _ 0.14 

\A DCS (Dj K°K+)/A CF (Df -► K°K+)\ 


(c) Black: la range predicted in Ref. m 


(d) Black: lcr range predicted in Ref. ED 


0.010 0.012 0.014 0.016 0.018 

B(d; -> k l k+) 

(e) 


FIG. 9: Blue: Our results for several observables probing doubly Cabibbo-suppressed amplitudes (see Eqs. ( |49| ) and ( |50[ )). The 
lines correspond to la (dashed), 2cr (solid) and 3<r (dashed-dotted) confidence intervals, respectively. The experimental error 
in R(D°) is obtained by Gaussian error propagation from Tab. [x] The results from other groups are shown in black. 

In case of R(D°) (a) and | ^cre^ol^^o^o) | (b) no errors are given in Refs. [52H55) . 
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Observable 

Measurement 

References 

SCS branching ratios 


B(D° - 

-y K+K~) 

(3.96 ± 0.08) • 10 -3 

m 


B(D° 

1 

+ 

t 

(1.402 ± 0.026) ■ 10 -3 

m 


B(D° 

-*• K S K S ) 

(0.17 ± 0.04) • 10 -3 

BU 


B(D° 

i 

O 

o 

(0.820 ± 0.035) ■ 10 -3 

m 


B(D+ 

-P 

o 

t 

(1.19 ± 0.06) • 10 -3 

Eel 


B(D+ 

->■ KsK+) 

(2.83 ± 0.16) • 10 -3 

m 

B(D+ - 

-+ K s n + 

)/B(Df -y 

K S K + ) (8.12 ± 0.28) • 10 -2 

EH 

B(D+ - 

o 

+ 

t 

)/B(D+ -y 

K S K+) (4.2 ± 1.4) ■ 10 -2 

m 

CF branching ratios 


B(D° 

-y K~tt + ) 

(3.88 ± 0.05) • 10 -2 

m 


B(D° 

K S n°) 

(1.19 ± 0.04) • 10 -2 

m 


B(D° 

-y k l k°) 

(1.00 ± 0.07) • 10 -2 

EH 


B{D+ 

— y K s ir+) 

(1.47 ± 0.07) • 10 -2 

EH 


B(D+ 

-»• K L TT + ) 

(1.46 ± 0.05) • 10 -2 

EH 


B(D+ 

-y K S K+) 

(1.50 ± 0.05) • 10 -2 

*13711381 

DCS branching ratios 

B{D° - 

+ 

* 

■f 

)/B(D° -y 

K~tt + ) 0.00349 ± 0.00004 

EH 


B(D+ 

-y K+ 77°) 

(1.83 ± 0.26) • 10 -4 

1351 



K + n~ 

strong phase difference 




Sk^ 

(6.45 ± 10.65)° 

*1111 


TABLE X: Input data for charm meson branching ratios and 
the strong phase difference 5kit used in our fit. Note that 
as we incorporate the correlations reported in Ref. [46J, for 
consistency we do not take into account the experimental up¬ 
dates of the following branching fractions: B(D° —> K~n + ) 


see Tabs. 
B(D+ -a 


XI 

XII 

STTZ 

and 

XIII 


^^ Note that B(D+ —y K + n°) and 
Ksn + ) are not part of the PDG fit, i.e., there are no 


correlation coefficients given for these decay modes. We there¬ 
fore have no correlation matrix for D+ decays. The value for 
B(D° —> I\ + n~)/B(D° —> K~tt + ) is taken from the Heavy 
Flavor Averaging Group (HFAG) in or der to take its corre¬ 


lation with Sktt into account, see Tab. 
*Our symmetrization of uncertainties. 


XIII 


*Our average. 


Appendix A: Input Data 

We give the input data used in the fits, including the 
correlation coefficients, in Tabs. [X] | X1V | For details on 
the input values for the form factors see Appendices |C 1| 
and 1C 21 


Appendix B: Mapping of the topological on the 
SU(3)p parameterization 


As discussed in Sec. |IIC| the topological flavor-flow pa¬ 
rameterization and the linear SU(3) f expansion can be 
mapped onto each other after the removal of redundan¬ 
cies in each parameterization. In Tabs. |XV| and |XVI| 
we give two numerical examples for the mapping. Both 
redundant parameters and redundant decay amplitudes 
have to be removed in order to obtain two corresponding 
11 x 11 regular coefficient matrices. Then, the mapping 



k>K+K~ 

£>£>o 

KJ7T ' 7T 

£> D o 

p? K S K s 

°£>° 

fO7T 0 7T 0 

toK 7r"*~ 

£> d o 

K k s ^ 
°£)0 

yoK+K- 

1.00 

0.38 

0.03 

0.09 

0.60 

0.21 

Bl>~ 

0.38 

1.00 

0.03 

0.09 

0.62 

0.22 

b k d s a Ks 

0.03 

0.03 

1.00 

0.01 

0.05 

0.03 

Y2 7T 0 7r° 

o D 0 

0.09 

0.09 

0.01 

1.00 

0.14 

0.05 

yoK~ ~k' 

D d o 

0.60 

0.62 

0.05 

0.14 

1.00 

0.35 

f2 K S^° 

°DO 

0.21 

0.22 

0.03 

0.05 

0.35 

1.00 


TABLE XI: Correlation coefficients for D° branching ratios 
[461 used in our fit. We abbreviate B{ = B(i -A /). Note 
that B(D° —y Klf°) is not part of the PDG fit and is used 
without correlations to the other modes. 



„K S K+ 
°D + 

k K s tt+ 
°D + 

k K + it° 
C D+ 

°d+ 

1.00 

0.75 

0.05 

n A 's 7r+ 

°D + 

0.75 

1.00 

0.06 

K K+n° 

0.05 

0.06 

1.00 


TABLE XII: Correlation coefficients for D + branching ratios 
[46] used in our fit, see the caption of Tab. XI for the notation 
used. Note that B(D + —> Klf + ) and B(D^ —> 7r°7r + ) are not 
part of the PDG fit and are used without correlations to the 
other modes. 


can be calculated by inverting one or the other coefficient 
matrix. We choose to omit the redundant amplitudes 

D° -A 7r°7T 0 , D° -A K~ 7T+ , D° -A I< + n~ , (Bl) 

D° -A K°tt° , D+ -A K°tt + , D + -A I< + ir° , (B2) 

using the sum rules presented in Appendix [D] We next 
calculate the redefined SU(3)jr matrix elements in terms 
of the remaining decay amplitudes by inverting the 
SU(3 )f coefficient matrix given in Tabs. I and V of 
Ref. [23] . The result for this inverse matrix is given in 
Tab. |XVlT| In order to illustrate how to read Tab.[XVH] 
we exemplify 

^27 = ^ A ( D+ /?07r+ ) + K ° K+ ) ' 

(B3) 

Inserting the expansions of A(D + -a K°n + ) and 
A(D f -A K°K + ) in terms of topological amplitudes into 
Eq. ( |B3| ) gives the desired expression of the SU(3 )_f ma¬ 
trix element in terms of the topological amplitudes. 
Note that in the matching we implicitly disregard higher 
order SU(3)i7’-breaking effects which are included in the 
approximate factorization formulas. Strictly speaking, 
these invalidate the linear SU(3 )f sum rules, see also 
Sec. |III A| However, this can be safely neglected as we 
only aim at a description of the data at linear 811(3)^ 
breaking here. 

While the exemplified topological-amplitude fit points 
respect the SU(3 )_f power counting, the SU(3 )_f break- 
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3ktt 

B(D u ^K + tt~) 

B(D°^K-tt+) 

Sktt 1.000 

0.404 

B(D°^>K-tt+) U ' 4U4 

1.000 


TABLE XIII: Correlation between B(D° ->• K + f~)/B(D° -»• 
K~ tt‘ ) and fevr El] used in our fit. 


F 0 DK ( 0) 

0.737 ± 0.005 

+ insHszi 

Fo n (0) 

0.638 ±0.012 


F 0 D ‘ K ( 0) 

(1 ± 5%) x F^( 0) 

EU [69] 


TABLE XIV: Numerical input for the form factors. The form 
factor F^ >aK (0) is varied flatly within the theory uncertainty 
m- Table adapted from ED- + Our average. 


ing matrix element s can nevertheless be quite large, like 
141 ~ 0.7 in Tab. |XV|in case of example point I. This 


shows that several small SU(3 )f breaking parameters 
of the topological-amplitude fit can add up to a larger 
SU(3 )f breaking matrix element of the group-theoretical 
approach. However, as demonstrated by the example 
point II in Tab. |XVI[ there are also solutions where both 
diagrammatic and group theoretic languages give SU(3 )f 
breaking < 50%. 


Appendix C: Approximate factorization formulas 

Below, we give the l/A^-leading expressions for the 
tree and annihilation diagrams. Corrections of higher 
order in the l/./V c -expansion are parameterized by St and 
8a introduced in Sec. |III A[ 


1. Factorization of tree amplitudes 


We use the following expressions for the l/iV c -leading 
contributions to the tree diagrams. SU(3 )f breaking in 
the l/iV 2 corrections is of higher order in our power 
counting and neglected, i.e., we use a flavor-universal 
correction parameter St- In our fit we vary 

0< |<$ T | < 0.15T fac , (Cl) 

—7r < arg(<$T) < 7r, (C2) 


with T fac defined in Eq. (24) and St = T — T fac , see 
Eq. (23). The l/AT c -leading, factorizable contributions 
to the SCS tree amplitudes are altogether given as: 


Tdo^ k +k- ~ k{^ 2 d Wr) K ( m ~ K ), (C3) 

T D°^n+ 7T- = fA™b - ml)F^(ml ), (C4) 


Topological parameter 

Value 

SU(3)f matrix element 

Value 

\Sa\/T^ 

0.14 


0.32 

\5t\/T fac 

0.15 

|4 5 | 

0.22 

arg(<5 A /T fac ) 

0.17 

141 

1.00 

arg (St/T^) 

-3.06 

|4| 

0.67 

\F 0 K ”(ml M )\ 

3.54 

|4| 

0.22 

F 0 DK ( 0) 

0.74 

i4h 

0.36 

EcP"(0) 

0.64 

l4 5l l 

0.39 

f o d ‘ k (0)/f o d «(0) 

0.95 

I4 52 l 

0.29 

arg (F^(m\ a) )) 

-1.74 

\Bl^\ 

0.18 

|C/T fac | 

1.10 

\Blr 2 \ 

0.07 

\E/T iac \ 

0.46 

l4 4l l 

0.13 

|Pb rea k/T fac | 

0.05 

arg(A^) 

1.44 

arg(C) 

2.47 

arg (4*) 

-2.53 

arg(c4) 

-1.50 

arg(A|) 

0.20 

arg (C«) 

-1.40 

arg(-Bi) 

-0.53 

arg(C 3 (1) ) 

0.00 

arg (4) 

-0.96 

arg (E) 

1.49 

arg(S 8 61 ) 

-2.11 

arg (44 

-0.65 

arg(4 51 ) 

-1.35 

arg(F^ 1) ) 

-0.93 

arg(4 52 ) 

-2.30 

arg(44 

-1.16 

arglPzr 1 ) 

2.56 

arg(Pbreak) 

0.18 

arg(4? 2 ) 

3.10 

\cP/c\ 

0.07 

arg (4 41 ) 

0.00 

\cP/C\ 

0.16 



\c^/c\ 

0.19 



\eP/e\ 

0.50 



\eP/E\ 

0.50 



\e£ ) /E\ 

0.05 



4‘° po 

0.50 



c/,C 

S x 

0.50 



5 X 

0.31 



X^^break 

°X 

0.07 



x 2 

0.27 




TABLE XV: Fit example point I and corresponding point in 
the SU(3 )f decomposition with linear SU(3 )f breaking. The 
values quoted for A) and B] in the last column are normalized 
to the largest SU(3 )f limit matrix element. The fit is only 
sensitive to St and 5 a in the combination (St — 5 A)/T ia,c = 
0.29e _3 02i . 


/ 7 -ifac _ Gp 1 f / 2 2 \ ttiDit ( 2 \ 

^Z)+—>7T+7r° — ~^ ; ~y^, ai ^ 7T ^ rnD ~ \ m TT ) ' 


(C5) 


G f 


T^ k+R o = - m 2 K )F 0 DK (m %), (C6) 


T 


tfac 


G f 


d+^tt+ko = -^ a ifAm Ds - m K )F 0 ‘ (m n ).( C7) 


The l/7V c -leading, factorizable contributions to the CF 
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Topological parameter 

Value 

SU(3)f matrix element 

Value 

\Sa\/T Uc 

0.15 

4271 

0.35 

|<H/T fac 

0.15 

41 s 1 

1.00 

arg (Sa/T^) 

1.27 

411 

0.19 

arg (5 t /T^) 

-1.88 

4?l 

0.36 

\F^(m 2 D J\ 

4.50 

411 

0.10 

F 0 DK ( 0) 

0.74 

4s Sl l 

0.15 

F?*(0) 

0.64 

4s 1Sl l 

0.49 

F o d ° k (0)/F o d ”(0) 

0.95 

4s 152 1 

0.06 

arg (F^(ml M )) 

-1.50 

427 5l l 

0.17 

|C/T fac | 

1.17 

427 52 l 

0.21 

|F/T fa0 | 

2.05 

4 2 2 7 4i I 

0.06 

|P br e a k/T faC | 

0.39 

arg(A^) 

0.94 

arg(C) 

2.23 

arg(Ag 5 ) 

-0.60 

arg(c) 1) ) 

-1.06 

arg(A|) 

1.47 

arg(C^ 1) ) 

1.48 

arg(Si) 

2.15 

arg(C^ 1) ) 

0.22 

arg(-Bf) 

0.55 

arg (E) 

2.57 

arg(S 8 61 ) 

-2.04 

arg (E^) 

0.23 

arg(B 8 1Sl ) 

-2.00 

arg (E^) 

1.72 

arg(B 8 152 ) 

-1.25 

arg (Ei 1) ) 

0.75 

arg (B 27 1 ) 

-2.84 

arg(Pbreak) 

1.97 

arg4a7 2 ) 

-1.53 

\c^/c\ 

0.22 

arg (B 27 1 ) 

0.23 

\cP/C\ 

0.38 



\C^/C\ 

0.12 



\E[ 1) /E\ 

0.06 



\E^/E\ 

0.09 



\eP/e\ 

0.31 



4 topo 

0.50 



ct,C 

S X 

0.50 



c f ,E 

5 X 

0.16 



Jf'^break 

°X 

0.50 



X 2 

0.12 




TABLE XVI: Fit example point II and corresponding point in 
the SU(3 )f decomposition with linear SU(3 )f breaking. Cf. 
Tab. [XV] for the notation. The fit is only sensitive to St and 
<Sa in the combination ( St — <Lt)/T fac = 0.30e -1 ' 88 *. 

tree amplitudes are given as: 

TdZ-+ K -*+ = ^| a i U( m D ~ m 2 K )Ff K (rnl ), (C8) 

t d+^k°„+ = ^= a i U(m 2 D - m 2 K )Ff K (ml). (C9) 

The l/iV c -leading, factorizable contributions to the 
DCS tree amplitudes are given as: 

fK(m 2 D - ml)F^ n (m 2 K ), (CIO) 


Td+^k+t o = - ml)Fo*(m 2 K ) , 

(Cll) 

T d+^k°k+ = -m 2 K )F°‘ K (m 2 K ). 

(C12) 


The matrix element of the vector current can be param¬ 
eterized by the vector and scalar form factors as [65] 


(P\ V»\D) =F°^ p (q 2 ) 


u 

Pd 


Pp 


m D ~ m \ 


+ Fr 


D^K 


( 9 2 )- 


m D — m 


K 


(C13) 


with the vector form factor F p ~ >p and the scalar form 
factor F P ~ J,K obeying [551 


(P\S\D)=F^ p (q 2 ) 


2 2 
m D — rrip 

m c — m p 


(C14) 


Here the same renormalization scheme and scale must be 
used for S and m c — m p . 

We calculate the form factors that appear in the tree 
amplitudes using the overall scaling factor appearing in 
the 2 -parameterization, i.e., a pole factor [mm 


f dk, m 2\_ F o >K (°) 

° [ p) 1-774/771^(23174’ 

(C15) 

F DsK (m 2 1 - F o‘ K (°) 

0 [ P> 1 -77l|/77l| ) ,(2400) ± ’ 

(C16) 

,... r."(«) 

0 p 1 — m 2 p /m 2 D , (24004 

(C17) 

with the scalar resonances [35] 


tof)» o ( 2317) ± = (2317.8 ± 0.6) MeV, 

(C18) 

77177 * (2400) ± = (2403 ± 40) MeV. 

(C19) 

The used input values for F PK ( 0), F PsK (0) and T/P^O) 
are given in Tab. XIV 

As we assume isospin symmetry in the topological- 


amplitude decomposition we ignore the smallish isospin 
breaking between charged and neutral masses of kaons 
and pions for consistency. We use the neutral masses in 
all amplitudes. However, in the phase space factors of 
the branching ratios we take the isospin mass splittings 
into account. 


2. Factorization of annihilation amplitudes 


We use the following expressions for the 1/TVc-leading 
contributions to the annihilation diagrams. As in the 
case of tree amplitudes we vary 6a of Eq. (25) as 


0 < \6 a \ < 0.15T fac , (C20) 

—7r < arg(<4) < 7r . (C21) 






















The l/iV c -leading, factorizable contributions [33] to the 
SCS annihilation amplitudes are given as: 


22 


a fac 

J± D+-+K°K+ 


= o, 


a fac 

^ 1 D+->if°7r+ 


Gf 

71 


aifD s F^ n {m 2 Ds ) ( m 2 K 



(C22) 

(C23) 


a fac 

A Dt- 


>K+ 


Gf 1 

r o — /— /^Q'ljDs 


V2 V2 


7 T'K'K 
^0 


(to 




(C24) 


The l/./V c -leading, factorizable contribution to the CF 
annihilation amplitude is given as: 

Ai Dt^KO K+ =°- ( C25 ) 

The l/iV c -leading, factorizable contributions to the DCS 
annihilation amplitudes are given as: 

a d+^k°*+ = f D F^{m 2 D ) (m 2 K - ml) , (C26) 

A%+^ K+n o = -^j=^=a 1 f L ,F^{m 2 D ) (m\ - m 2 ) . 

(C27) 


Note that the l/IV c -leading SCS annihilation amplitude 
^ d + k°k+ an< ^ the l/iV c -leading CF annihilation am¬ 
plitude can be neglected due to isospin sym¬ 

metry m■ However, the corresponding 1 /N% corrections 
are of course taken into account (as for the others) and 


specified in Tab. VI 


Constraints on \Fq ^ ( 777 z? ( ) I can be taken from r de¬ 
cays. In order to accommodate the measurements of 
r — > KgTT~v T from Belle [3H] we vary the form factor 
in the interval 


l<\F 0 Kn (rn 2 D J\<^, (C28) 

-7r<arg(F 0 ^(m|, (s) )) <tt, (C29) 

setting 

F?”(Tn 2 D.)=F?”(m 2 D ). (C30) 


Appendix D: Diagrammatic representation of sum 
rules 


In Tabs. XVIII XXIII we give the diagrammatic rep¬ 
resentation of the six Grossman-Robinson SU(3 )f sum 
rules which hold to linear order in SU(3 )f breaking [29] , 
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Decay amplitude 

T 

c 

E 

-fbreak 

+ -^A(D°-> *+*-) 

1 

X 

Ik 

0 

-T» s=<g 


+A(D° 7r°7r°) 

0 


+ 73* 


~A(D + 7T°7r + ) 



0 

0 


+7 3 x(r^i 


TABLE XVIII: Diagrammatic representation of sum rule I, AjA (D° n + n ) + A(D° —> n°n°) — A(D + —> n°n + ) = 0. 


Decay amplitude 

T 

T a) 

c 

cf> 

E 


^A(D°^K~n+) 


1 , 

0 

0 



75 tri 

- x O 

+3*E=<^ 


+A(D° ->• K°n°) 

0 

0 

75* f) 

* X (T^ 

7!* 

tL C<S 

-73 A(D+^K°n+) 

, 


1 

i yft _^ 

0 

0 

~' /I ( 14 ) 

7S>< 0 

^ rn) 

~^ x r ^n 




TABLE XIX: Diagrammatic representation of sum rule II, -^A(D° —> K 7r + ) + A(D° —> K°tv°) — ^A(D + —> K°n + ) = 0. 


Decay amplitude 

T 

7^(1) 

1 2 

A 

4 1J 

c 


E 

rj 2 

A{D° -»• A+tt”) 

1 *(T J8 ^ 

u . u 

lxp^j 
u . u 

0 

0 

0 

0 

lx 

lx (D<2 

V2A{D° ->■ K 0 tt°) 

0 

0 

0 

0 

lx /p-^ 

lx n ■ 

U—i—U 

-lx ocg 

-lx DCg 

-A{D+ ->■ if°7T + ) 

0 

0 

'* ><A 

'* (X$ 

— lx p-<^ 

— lx n ■ ^2? 

U—.—U 

0 

0 

-V2A(D+ -»• A+tt 0 ) 

-lxp^ 
u U 

-lxp^ 

u , u 

lx 3>< 

lx (X£ 

0 

0 

0 

0 


TABLE XX: Diagrammatic representation of sum rule III, A(D° — y K + tt ) + \/2A(D° — > K°tv°) — A(D + — ¥ K°n + ) — 
V2A(D+ -> K+n°) = 0. 
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7 T + 7 r 


TABLE XXI: Diagrammatic representation of sum rule IV, A(D° K K + ) — A(D° 
K+ tt-)=0. 


) - A(D° -»• K~ 7T+) - A(D° ->• 
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T - ' 
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